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Charge carrier solvation and large polaron formation on a polymer chain revealed in 

model ab initio computations 

M. L. Mayo^ and Yu. N. Gartstein^ 

'Department of Physics, The University of Texas at Dallas, 
P. O. Box 830688, F023, Richardson, Texas 75083, USA 

When an excess charge carrier is added to a semiconducting polymer chain, it is well known that 
the carrier may self-trap into a polaronic state accompanied by a bond length adjustment pattern. 
A different mechanism of self-localization is the solvation of charge carriers expected to take place 
when the polymer chain is immersed in polar media such as common solvents. We use state-of-the- 
art ab initio computations in conjunction with the Polarizable Continuum Model to unequivocally 
demonstrate solvation-induced self-consistent charge localization into large-radius one-dimensional 
(ID) polarons on long CivH2 carbon chains with the polyynic structure. Within the framework used, 
the solvation results in a much more pronounced charge localization. We believe this mechanism of 
polaron formation to be of relevance for various ID semiconductors in polar environments. 

PACS numbers: 31.15.A-, 31.15.at, 71. 15. Mb, 71.38.Fp, 81.07.Nb 
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A variety of one-dimensional (ID) wide-band semicon- 
ductor (SC) structures: conjugated polymers, nanowires, 
and nanotubes, attract a great deal of attention. The 
nature and properties of excess charge carriers on these 
structures are fundamental for various applications. For 
our discussion, one distinguishes between nearly free 
band states of excess carriers and self-localized polaronic 
[l| states. Self-trapping and the formation of ID polarons 
can occur owing to the (strong) interaction of the elec- 
tronic subsystem with another subsystem such as, e.g., 
displacements of the underlying atomic lattice, which is 
the mechanism that has been extensively explored for ID 
electron-phonon systems. Here we are concerned with a 
different, much less studied implementation of the pola- 
ronic effect that should take place for ID SCs immersed 
in 3D polar media, the situation particularly common 
for applications involving fundamental redox processes 
in polar solvents. In this case of what could be called 
charge carrier solvation, the long-range Coulomb inter- 
action is expected 0, y, lj| to lead to the formation of a 
localized electronic state on a SC structure surrounded 
by a self-consistent pattern of the sluggish (orientational) 
polarization of the solvent. The main goal that we pursue 
and achieve in this paper is to substantiate the validity 
of this picture at the level of ab initio computations, that 
is, with realistic many-electron effects taken into account. 
To our knowledge, this is a first ab initio demonstration 
of the solvation-induced self-consistent localization of the 
excess charge resulting in large-radius ID polarons. 

The type of systems we consider is an interesting model 
realization in the general context of quanturnparticles in- 
teracting with dissipative environments [5|, |6| . The con- 
finement of the electron motion to ID and its structural 
separation from the 3D polarizable medium clearly dis- 
tinguishes such systems from the well-known 3D cases of 
polarons in polar SCs [71, Is], Isf and solvated electrons in 
polar liquids [9|, llO|. A combination of several factors 
here creates conditions that may be especially favorable 



for the formation of ID continuum adiabatic polarons; 
these are small effective band masses of the carriers (e.g., 
~ 0.1 of free electron mass), enhanced Coulomb effects in 
ID [111 , and the slow response of the solvent polarization. 
The self-trapping of an excess charge carrier (an elec- 
tron or a hole) is easily described in a single-particle pic- 
ture within a standard ID continuum adiabatic frame- 
work. It corresponds to the localized ground-state solu- 
tion iIj{x) of the non-linear Schrodinger equation: 



2m dx^ 



dx'V{x - x') l-ipix')]^ ipix) = E%l;{x), 



where x is the coordinate along the structure aixs, m the 
intrinsic effective mass of the carrier, and V{x) the effec- 
tive self-interaction mediated by another subsystem. In 
the case of a short-range electron-phonon mediation, the 
self-interaction can be taken local: V{x) = gS{x) lead- 
ing then to the exact result ipi^) o(. sech.(gmx/2h^) for a 
continuum ID polaron widely known after the pioneer- 
ing contributions of Rashba [12] and Holstein [l3|. In 
the case of the long-range polarization interaction, the 
effective V{x) behaves as l/\x\ at large distances, while 
the short-range behavior depends on specific system de- 
tails such as, e.g., the actual transverse charge density 
distribution and the geometry of the dielectric screen- 
ing. Numerical studies of this case [3l, |4| indicate that 
the binding energy of the resulting polarons could reach 
a substantial fraction, roughly one-third, of the binding 
energy of the corresponding Wannier-Mott excitons (pri- 
mary optical excitations in many ID SCs). 

While illuminating the basic physics of the effect, the 
single-particle description has the drawback of not ex- 
plicitly including valence band electrons, whose reorga- 
nization may substantially affect the outcomes. Such a 
reorganization can be important even in non- interac ting 



electron models - a nice illustration was given in Ref. 14 1 
that analyzed how polarons of a two-band Peierls dielec- 
tric model evolve into single-particle Holstein polarons 



in the limit of the "frozen valence band" approximation. 
With realistic Coulomb interactions in place, the role of 
valence electrons in the formation of the relevant self- 
consistent potentials would only increase. As applied to 
electron-phonon polarons in conjugated polymers (CPs), 
various ab initio computational schemes have been found 
useful to clarify the role of electron-electron interactions 
on such polarons and their optical properties (we refer 
the reader to a recent account [15| for a discussion and 
multiple references). In what follows we use the Gaus- 
sian 03-implemented state-of-the-art DFT ab initio cal- 
culations in conjunction with the Polarizable Continuum 
Model (PCM) [3, [13] to study the polaron formation in 
conjugated linear systems as mediated by the polarization 
interaction with a surrounding polar solvent. 

Since the computational demand increases substan- 
tially in DFT-PCM calculations, we have selected struc- 
turally simple systems of hydrogen-terminated polyynic 
linear carbon chains CArH2 (even N between 20 and 100) 
for our demonstration. A polyynic structure exhibits in 
its ground state an alternating pattern of triple C = C 
and single C — C bonds and is semiconducting. It should 
be noted that, while playing the role of a prototypical 
example in our study, polyynic chains continue to be the 
subject of much attention in their own right [18Lll9l.l2Cll|. 

The results presented in this paper have been derived 
by employing DFT computations with the B3LYP hybrid 
exchange-correlation functional known to include local 
and non-local correlation effects. The 6-311-1— |-G(d,p) all 
electron basis set was chosen due to its richness and abil- 
ity to represent charged systems. We note that our "in 
vacuum" calculations have been verified to compare well 
against previously published ab initio results on neutral 
[l8| and charged Jl9] polyynic systems both in terms of 
energetics and optimized bond lengths. The DFT-PCM 
( "in solvent" ) calculations used water as a polar solvent 
with its default parameters in Gaussian 03. 

A pertinent quantity for us is the equilibrium spatial 
distribution of the excess charge over the ID SC, which 
is naturally related to atomic charges in outputs of ab 
initio calculations. Importantly, these charges represent 
the effect due to all electrons in the system. It is well 
known [i7| that calculations of the atomic charges from 
many-electron wave functions are dependent on the ba- 
sis set transformations and sometimes lead to spurious 
results. Two procedures, due to MuUiken and Lowdin, 
are used particularly frequently in such calculations [l7| . 
We exercised both and established that, while they may 
lead to different raw atomic charges for individual cases, 
the results become remarkably close when calculated for 
unit-cell-based charge differences between charged and 
neutral chains. The unit cells are defined here to consist 
of pairs of neighboring carbons, with the chain end cells 
including the hydrogen atomic charges. It is these stable 
charge differences - that is, excess charge - per cell that 
are displayed in Figs. [1] [2] and SJa) using Lowdin num- 
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FIG. 



1: Excess charge density spatial distribution on charged 
CatH^ chains with iV = 20, 40, 60, 80 and 100 - top to bot- 
tom curves, respectively, in their central parts. All chains 
feature the same pattern of fixed bond lengths corresponding 
to RGl geometry described in text. Panel (a) displays re- 
sults for chains in vacuum, panel (b) for chains in the solvent 
environment. 



bers. We studied both anionic (an extra electron) and 
cationic (an extra hole) chains. Their self-localization 
behavior has been found very similar. Accordingly we 
restrict our discussion here to the anionic cases, with the 
excess charge in the figures measured in units of (— e), 
where e is the magnitude of the fundamental charge. 

To be able to clearly separate the effect of the solvation 
on excess charge localization from the effect of atomic 
displacements, we start by exploring systems with rigid 
geometries, that is, with fixed prescribed atomic posi- 
tions. Our first geometry, referred to as RGl, features 
the following fixed lengths: 1.235 A for the triple bonds, 
1.326 A for the single, and 1.063 A for the C - H bonds. 
These lengths have been retrieved from our complete op- 
timization of the neutral N = 80 chain in vacuum and 
are therefore expected to be relatively close to optimal 
values. We now use this geometry to perform DFT calcu- 
lations for a set of C7vll2 and CatHj" samples with various 
N to find the excess charge density distribution both in 
vacuum and in the solvent environment. The derived re- 
sults are compared in Fig. [1] whose (a) and (b) panels 
display clearly distinct trends. Barring the end effects, 
vacuum results in panel (a) are reminiscent of the "par- 
ticle in a box" behavior: as the chain length increases, 
the excess charge is distributed more and more uniformly 
over the whole chain. In a sharp contrast, for chains in 
the solvent, panel (b), the excess charge exhibits much 
more localized distributions around central parts of the 
chains. As the chain length increases, those distributions 
show a tendency towards convergence, albeit not com- 
pletely achieved within the range of lengths studied. 

To further the understanding of the observed behav- 
ior, we now attempt to decrease the spatial extent of 
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FIG. 2: Excess charge density spatial distribution on charged 
CatH^ chains, (a) N=80. The evolution of the charge distri- 
bution upon stretching the single bonds by 0.1 A successive 
increments from 1.326 A to 1.726 A - bottom to top in a 
group of 5 clearly distinct curves shown by different connected 
symbols and denoted "in solvent". The same - but uncon- 
nected - data symbols are used to show the results for chains 
in vacuum, those distributions largely overlap, (b) Compar- 
ison of distributions for chains in the solvent with A*' = 40, 
60, 80 (different data symbols) and 100 (solid line) for the 
geometry RG2 featuring the longest single bonds. 



charge localization. In single-particle models, the size 
of the polaron is determined by the balance of the gain 
in the potential energy of the localized carrier and the 
loss in its kinetic energy. Decreasing the latter leads to 
a diminished delocalization propensity and is expected 
to shorten the polaron. This should be achievable with 
an increased intrinsic effective mass, or with a narrower 
electron bandwidth. To imitate the effect, we chose to ar- 
tificially stretch all single bonds in our system, thereby 
making C = C "dimers" more and more separated from 
each other. Figure ^a,) illustrates a dramatic difference 
in the response of the excess charge distribution to suc- 
cessive single-bond stretchings for chains in vacuum and 
for chains in solvent. Whereas charge density in vac- 
uum is hardly affected in this process, preserving nearly 
uniform distributions over the whole chain, for chains in 
solvent, each successive stretching step results in smaller 
and smaller localization regions. We refer to the fixed ge- 
ometry with the longest-stretched single bonds of 1.726 
A as RG2. As the localization size in this geometry ap- 
pears to be well within the chain length, we compare the 
excess charge distributions for several chain lengths in 
Fig. mb). The distributions for N = 60, 80, and 100 
practically coincide confirming that they correspond to 
a fully converged self-consistent pattern of excess charge 
localization achieved entirely due to the interaction with 
the surrounding solvent. 

A comparison of the corresponding molecular orbital 
(MO) energy level structures is available in Fig. [31 One 
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FIG. 3: Molecular orbital energy levels for neutral (0) and 
charged (-1) C100H2 chains with rigid RGl (a-d) and RG2 (e- 
h) geometries both in vacuum and in the solvent. The arrows 
indicate positions of HOMO levels. 



immediately notices a much wider (but inconsequential 
for the effect we consider) band gap in the RG2 geometry. 
The narrowing of the bandwidth is, however, relevant 
and easily detected via a much larger density of states 
(panels (e) vs (a) and (g) vs (c)). Disregarding "trivial" 
(the electrostatic potential with respect to infinity) over- 
all shifts of the levels for the charged chains in vacuum, 
the important difference with the charged chains in the 
solvent is that the latter case exhibits local highest occu- 
pied (HOMO) levels clearly separated from other states: 
compare panel (d) to (b) and, especially, panel (h) to (f). 
This is the behavior expected from single-particle polaron 
models. One, however, does not observe distinctly sep- 
arate unoccupied local levels appearing in such models 
with long-range Coulomb potentials [3J. Also, a reorga- 
nization of valence band electron states is apparent in the 
band gap upon the formation of the polaron, especially 
well seen in geometry RG2, panel (h). This is suggestive 
of a more complex electronic structure of the polaron as 
affected by many-electron interactions. 

To better benchmark the localization due to the sol- 
vation, we now turn to bond-length optimization results 
derived for the A^ = 40 system. Figure|3Ua) pictures mod- 
ifications of the excess charge spatial distribution that 
take place upon changes from the rigid RGl geometry 
to fully optimized bond length patterns both in vacuum 
and in the solvent. It is evident that for the charge dis- 
tribution, the effect of the lattice relaxation is much less 
significant than the effect of the solvation. (The same 
conclusion has been drawn in our studies of iV = 80 
chains.) Figure |4{b) shows how local HOMO and intra- 
gap valence levels are affected by the pure solvation and 
by the combined action of the solvation and lattice dis- 
placements. Interestingly, here the effects due to each of 
the mechanisms appear similar in magnitude despite a 
much more pronounced effect of the solvation on charge 
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FIG. 4: Comparison of the solvation and bond-length adjust- 
ment effects for C40H2. (a) Excess charge distributions in 
vaccum (a) and in the solvent (D) for RGl geometry (dashed 
lines) and for the fully optimized bond length patterns (solid 
lines), (b) The corresponding MO energy levels in the sol- 
vent for: (bl) Optimized neutral chain, (b2) Charged chain 
in RGl geometry, (b3) Optimized charged chain. The arrows 
show HOMO levels. 



density localization. This is again likely indicative of the 
intricacies of many-electron interactions and the screen- 
ing effects. A more complete understanding of similari- 
ties and differences of these mechanisms calls for further 
studies at various computational levels. The picture of 
Fig. mb) suggests that comparative spectroscopy of elec- 
tronic transitions to and from local HOMO levels should 
be a helpful experimental tool. 

The spatial extent of demonstrated solvation-induced 
ID polarons is seen to be on the order of tens of 
Angstroms. We, however, do not feel to be able to reli- 
ably evaluate the polaron binding energy except to esti- 
mate it at ^ 0.1 eV based on comparisons of LUMO levels 
with the computed electron affinities. Our consideration 
has been limited to chains with idealized straight struc- 
tures; various defects including spatial conformations of 
the chains can lead to the formation of bound polarons 
with possibly larger binding energies. 

While providing a clear conceptual demonstration, the 
present ab initio framework may overestimate the effects 
of the solvation on self-localization as it does not take a 
proper account of the frequency dependence of the sol- 
vent dielectric function [3|, |j, [Zl, |9|| . This issue could possi- 
bly be addressed within time-dependent DFT-PCM com- 
putational schemes 



21|. Another current development. 



on the scaling description of solvent effects 22] , may also 
help in studies of large polaron formation. Given the gen- 
eral character of the solvation-induced self-localization 
with the subsequent decrease of charge carrier mobility 



[2, la l23|, as well as its expected significance for fun- 
damental redox processes on various ID SC structures, 
further studies of this effect seem quite interesting and 
warranted. 
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